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Method and System for Free Space Optical Tomography of Diffuse Media 

CROSS REFERENCE TO RELATED APPLICATIONS 

This application claims the benefit under 35 U.S.C. §1 19(e) of U.S. Provisional 
5 Application No. 60/445,016 filed February 5, 2003, which application is incorporated herein by 
reference in its entirety. 

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH 
Not Applicable. 

10 

FIELD OF THE INVENTION 

This invention relates generally to tomography and, more particularly, to optical 
tomography. 

15 BACKGROUND OF THE INVENTION 

As is known in the art, tomography is often used in construction of detailed images of 
internal structures of objects. Tomography relies upon a selected form of energy being directed 
toward and passing through the object at more than one angle. The energy from the various 
angles is collected and processed to provide a tomographic image. The received signals are 
20 typically less intense (i.e., darker) where the object is thicker or more dense, and more intense 
(i.e., brighter) where the object is thinner or less dense. 

As is known, a signal received by a single energy sensor (i.e., at one angle) does not 
contain sufficient information to generate either a two-dimensional or a three-dimensional 
representation of internal structures of the object. As is also known, signals received by energy 
25 sensors arranged in a plane or volume provide sufficient information to generate a three- 
dimensional representation of internal structures of the object. 

Tomography is used in a variety of systems with a variety of types of transmitted and 
- received energy. For example, in x-ray Computed Axial Tomography (CAT), x-ray energy is 
projected through an object, typically at a variety of angles, and a variety of x-ray receivers, at a 
30 corresponding variety of angles, are used to receive the x-ray energy. A computer is used to 
generate an image of internal structures of the object in three dimensions from signals received 
by the variety of x-ray receivers. It should be recognized that x-rays tend to pass through the 
object in straight lines with relatively little attenuation. 
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One type of x-ray CAT system is used for medical imaging and the object through which 
the x-rays are projected is a person. However, the x-ray CAT system can be used to image 
internal structures of other objects, for example, luggage at an airport. 

Some forms of optical tomography are known, which use one or more wavelengths of 
5 visible or invisible light rather than x-rays. However, unlike x-ray tomography, for which x-rays 
tend to pass through an object in a straight line with relatively little attenuation, light tends to be 
absorbed and to scatter when passing though an object Therefore, light does not travel in 
straight lines when passing through the object. Light also tends to attenuate and to scatter more 
when passing though a relatively thick object having a relatively non-homogeneous medium, 

1 0 than it tends to attenuate and to scatter when passing through a relatively thin object having a 
relatively homogeneous medium. 

Diffuse Optical Tomography (DOT) and Fluorescence Molecular Tomography (FMT) 
are known optical imaging techniques that allow optical tomography imaging of internal 
structure of body parts of animals and humans. DOT is an effective imaging technique capable 

1 5 of imaging hemoglobin concentration and oxygenation. DOT can increase the specificity in 

detecting disease when used in combination with more established examinations (for example for 
cancer or arthritis detection and characterization). In addition, DOT is used to study brain 
activation and to study exercising muscle. 

FMT uses fluorochromes, which absorb light propagating inside of an object and emit 

20 light at a longer wavelength (lower energy) than the absorbed light inside of the object, allowing 
non-invasive in-vivo investigation of functional and molecular signatures in whole tissues of 
animals and humans. FMT enables molecular imaging, i.e., it can probe molecular abnormalities 
that are the basis of a disease, rather than imaging the end-anatomical effects of the molecular 
abnormalities as with conventional imaging approaches. Specific imaging of molecular targets 

25 provides earlier detection and characterization of a disease, as well as earlier and direct 

molecular assessment of treatment efficacy. FMT technology can also transfer typical in-vitro 
fluorescence assays, such as gene-expression profiling, elucidating cellular pathways or sensing 
small molecule protein interaction to in-vivo non-invasive imaging applications of large tissues. 
Some conventional optical tomography systems use near infrared (near-IR or NIR) light, 

30 instead of light in the visible spectrum when passing through animal tissues, since NIR tends to 
attenuate less than visible light. The use of NIR light instead of light in the visible spectrum 
provides the ability to image deeper tissues, i.e., thicker tissues, or with higher sensitivity than in 
the visible light region. The development of highly efficient fluorescent probes, i.e., 
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appropriately engineered fluorochromes with high molecular specificity emitting in the NIR, has 
also enabled FMT imaging of deeper tissues. 

Mathematical modeling of light propagation in animal tissue and technological 
advancements in light sources (photon sources) and light sensors (photon receivers or photo 
5 detectors) has made optical tomography possible using diffuse light Diffuse Optical 

Tomography (DOT) uses multiple projections and de-convolves the scattering effect of tissue. 

Conventional DOT and FMT systems include a light source (such as a diode laser and 
appropriate driver), and an optical switch, which provides light to a group of optical pathways, 
for example, optical fibers. The optical switch directs the light source to selected ones of the 
10 optical fibers, one at a time, in a sequence. The optical fibers are in direct contact with a diffuse 
medium to be imaged. Using the optical fibers, the single laser source is directed to selected 
points on the surface of the diffuse medium. Light is collected with the use of fiber bundles, 
placed at multiple points, also in direct contact with the surface of the diffuse medium, and the 
light is directed though the fiber bundles from the diffuse medium to appropriate light sensors. 
15 A computer performs tomographic data transformations to provide images for display and 
storage. 

Existing DOT and FMT systems employ light sources and light sensors in direct contact 
with the object to be imaged, providing direct contact systems. For a description of an 
exemplary optical tomography system, see D. J. Hawrysz and E. M. Sevick-Muraca, 
20 "Developments Toward Diagnostic Breast Cancer Imaging Using Near-Infrared Optical 
Measurements and Fluorescent Contrast Agents," Neoplasia, vol. 2, pp. 388-417, 2000. 

The contact light sensors each receive light essentially from a single respective point on 
> the surface of the object. Direct contact systems tend to reduce system versatility, limiting the 
shapes, sizes, and geometries of objects that can be tomograhically imaged with any particular 
25 DOT or FMT system. Direct contact systems, when used to image body parts of a patient, also 
tend to limit patient comfort. 

Some optical tomography systems can only be used to image objects having a particular 
shape. For example, Wake et al., U.S. patent number 6,21 1,512, describes an optical 
tomography system for use only with objects having a cylindrical geometry. 

30 

SUMMARY OF THE INVENTION 

In accordance with the present invention, a system and a method to generate an image by 
optical tomography of internal structures of an object, such as a body part, includes a surface 
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capture system to capture a surface model (e.g., a mathematical description) of the surface of the 
object and to provide the surface model to one or more optical models. The optical models can 
be used to account for propagation of the light toward light sensors, which can be spaced apart 
from the object. 

5 In accordance with the present invention, a system for optical tomography of an object 

includes a data collection system to collect light having traveled through the object and in free 
space about the object and to generate one or more signals in response to the light. The system 
also includes a surface-capture system to receive surface light reflected from a surface of the 
object and to generate a surface model (e.g., a three-dimensional mathematical description) of at 

10 least a portion of the surface of the object. The system further includes a model processor 

coupled to the surface-capture system to provide one or more optical models associated with the 
surface of the object. A solution processor coupled to the data collection system and the model 
processor provides a tomographic image of the object in response to the one or more signals and 
to the one or more optical models. In one particular embodiment, the data collection system 

15 includes one or more light sources for emitting light directed to propagate through the object and 
one or more light sensors disposed to receive portions of the light which have passed through the 
object, wherein at least one of the one or more light sources and the one or more light sensors are 
disposed apart from the object. 

In accordance with a still further aspect of the present invention, a method of imaging a 

20 diffuse object with optical tomography includes propagating light toward the object, capturing 
light, which travels through the object, and generating one or more signals in response to the 
light. The method also includes capturing surface light, which reflects off a surface of the object, 
generating a surface model (e.g., a three-dimensional mathematical description) of at least a 
portion of the surface of the object which reflected light, and generating one or more optical 

25 models associated with the surface of the object. The method further includes generating a 
tomographic image of the object in response to the one or more signals and to the one or more 
optical models. In one particular embodiment, the method also includes emitting light with one 
or more light sources, directing the light toward and through the object, and receiving the light 
with one or more light sensors. 

30 With this particular technique, a method of generating a tomographic image of an object 

using light energy is provided by generating a mathematical description of at least a portion of 
the surface of the object and using the mathematical description to generate optical models. At 
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least one of the one or more light sources and/or the light sensors can be spaced apart from the 
object. 

BRIEF DESCRIPTION OF THE DRAWINGS 
5 The foregoing features of the invention, as well as the invention itself may be more fully 

understood from the following detailed description of the drawings, in which: 

FIG. 1 is a block diagram of an exemplary optical tomography system; 

FIG. 2 is a block diagram of a portion of an exemplary optical tomography system, which 
includes a light sensor and a plurality of light sources; 
10 FIG. 3 is a block diagram of a portion of another exemplary optical tomography system, 

which includes a plurality of light sensors and a plurality of light sources; 

FIG. 4 is a flow chart showing an exemplary method for imaging an object with optical 
tomography; 

FIG. 4A is a flow chart showing further details of the method of FIG. 4; 
15 FIG. 5 is a flow chart showing another exemplary method for imaging an object with 

optical tomography; 

FIG. 5A is a flow chart showing further details of the method of FIG. 5; 

FIG. 6 is a pictorial used to describe an in-focus imaging system used in the optical 
tomography system of FIG. 1; 
20 FIG. 6A is another pictorial used to describe the in-focus imaging system used in the 

optical tomography system of FIG. 1; 

FIG. 6B is a pictorial used to describe the in-focus imaging system having a small 
'aperture and having a large aperture used in the optical tomography system of FIG. 1; 

FIG. 6C is a pictorial used to describe an isotropic light source used in the optical 
25 tomography system of FIG. 1 ; 

FIG. 7 is a pictorial used to describe an out-of-focus surface contribution received by a 
light sensor used in the optical tomography system of FIG. 1 ; 

FIG. 7A is a pictorial used to further describe an out-of-focus surface contribution 
received by a light sensor used in the optical tomography system of FIG. 1 ; 
30 FIG. 7B is a pictorial used to describe an out-of-focus a solid angle contribution received 

by a light sensor used in the optical tomography system of FIG. 1; 

FIG. 7C is another pictorial used to describe an out-of-focus solid angle contribution 
received by a light sensor used in the optical tomography system of FIG. 1; 

5 
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FIG. 7D is a pictorial used to describe an out-of-focus full surface contribution received 
by a light sensor used in the optical tomography system of FIG. 1; and 

FIG. 7E is another pictorial used to describe an out-of-focus contribution received by a 
light sensor used in the optical tomography system of FIG. 1. 

5 

DETAILED DESCRIPTION OF THE INVENTION 

Referring now to FIG. 1, an optical tomography system 10 includes a data collection 
system 12 having one or more light sources 14 spaced from an object 18 under test. Each of the 
one or more light sources 14 projects light 16 toward the object 18. Portions (not shown) of the 

10 light 16 which pass through the object 18 are received by one or more light sensors 20 which are 
disposed proximate, but spaced apart from, the object 18. It should be appreciated that sensors 
20 are disposed about the object 20 such that the sensors 20 can receive light, which propagates 
through the object 18. In one particular exemplary embodiment, the light sensors 20 are 
disposed to be approximately one centimeter apart from the object 18, proximate a side of the 

15 object 18 which is substantially opposite from the side of the object upon which the light 16 is 
directed. However, in other embodiments, the one or more light sensors 20 can be disposed 
more that one centimeter or less than one centimeter apart from object 18, and can be disposed 
apart from any side of the object. In another exemplary embodiment, the light sensors are 
disposed to be between five and twenty-five centimeters from the surface of the object. The 

20 separation between the surface of the object and the light sensors is selected in accordance with a 
variety of factors, including, but not limited to, a distance that can achieve a proper focal depth, 
while maximizing light collection capacity. 

The light sensors 20 receive the light, which passes through the object 18. Necessarily, 
since the one or more light sensors 14 are spaced from the object, the light propagates in free 

25 space prior to reaching the sensors 14. The one or more light sensors 20 provide corresponding 
light sensor signals (e.g. in the form of electrical or other types of signals) to a data processor 24. 
The data processor 24 digitizes, formats, and combines the digitized and formatted light sensor 
signals into vectors for subsequent processing, the functions of which are described in more 
detail in conjunction with FIGS. 4A and 5A. 

30 In some embodiments, the light sensors 20 are also adapted to receive fluorescent light 

generated by fluorochromes internal to the object 18, for example from fluorescent probes 
injected into the object 18 which tend to coalesce in particular structures or molecules within the 
object 18. 
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The data collection system 12 is coupled to a solution processor 40 and, in the exemplary 
embodiment of FIG. 1, the data collection system 12 provides measured optical data to the 
solution processor 40 through the data processor 24. The solution processor 40 provides a 
solution to an "image problem" described more fully below, which provides image data 
5 corresponding to internal structures in the object 1 8. 

The optical tomography system 10 also includes a surface capture system 28 having one 
or more surface light sources 30 spaced from the object 16. The one or more surface light 
sources 30 project surface light 32 at the object 1 8. Portions of the surface light 32 reflect from 
the surface of the object 1 8 and are received by one or more surface light sensors 34. The one or 

10 more surface light sensors 34 receive the surface light reflected from the object 18 and provide 
corresponding surface light sensor signals (e.g. in the form of electrical or other types of signals) 
to a surface capture processor 36. In response to the surface light sensor signals, the surface 
capture processor 36 generates a model (e.g., a mathematical description) of at least a portion of 
the surface of the object 18 from which the surface light reflects. As used herein, the term 

15 "mathematical description" can refer to an algorithmic description formed by equations, a 

numerical description formed by numbers, or both, and is also referred to herein as a model of 
the surface of the object. One surface capture system 28 generates the surface light 32 in a 
predetermined spatial pattern which is received by the one or more surface light sensors 34 
provided as one or more three-dimensional cameras. 

20 It will be recognized, however, that there are a variety of surface capture systems can 

generate a model of the surface of the object 18. 

The surface capture system 28 is coupled to a model processor 39. The model processor 
39 generates one or more optical models. The optical models are described more fully below. 
However, let it suffice here to say that, where the object 18 is diffuse to the propagation of the 

25 light 1 8 through the object, a first optical model generated by the model processor 39 can model 
the light 16 in the object 18. Furthermore, the first optical model can assume that the object 16 is 
not only diffuse but also that the object 16 is homogeneous vis-&-vis propagation of the light 16 
within the object. In other words, the first optical model can assume that the object 16 has no 
internal structures. Other optical models, referred to herein collectively as a second optical 

30 model, which can be generated by the model processor 39 include, but are not limited to, a model 
of the light 16 as it passes from inside of the object through the surface of the object, a model of 
the light 16 as it propagates through free space toward the one or more light sensors 20, and a 
model of optical characteristics of each of the one or more light sensors 16. 
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In the exemplary embodiment of FIG. 1, the model processor 39 provides the one or 
more optical models described above to the solution processor 40. Therefore, the solution 
processor receives the one or more light sensor signals (e.g., electrical signals) from the data 
processor and the one or more optical models from the model processor. 
5 The solution processor 40 is described more fully in conjunction with FIG. 4A. 

However, it should be understood that the one or more electrical signals provided by the data 
processor can correspond to measured optical data associated with the light 16 which has 
propagated through the object 18 and through free space before being collected by the one or 
more light sensors 20. It should be further understood that the one or more optical models 

10 provided by the model processor 39 can correspond to a theoretically derived "expected" 

response of the one or more light sensors 20 assuming that the object 1 8 is both internally diffuse 
and also homogeneous, e.g., having no internal structures. Where the object 18 does in fact have 
internal structures and is not internally homogeneous (for example, as is generally the case of a 
human body part), the solution processor 40 is presented with an "image problem" of the form: 

15 measurements = (theoretical predictions) x (unknown distribution), where the measurements are 
provided by the data collection system 12 and the theoretical predictions are provided by the 
surface capture system 28 in combination with the model processor 39. The solution processor 
40 can solve for the unknown distribution in order to establish physical positions and 
characteristics of the internal structures in the object 18. 

20 The solution processor 40 provides an output to an image processor 41 , which in turn 

provides data to a display 42. The display 42 can provide tomographic images of the internal 
structures of the object 18. 

As indicated by dashed line 38 in FIG. 1, in alternate embodiments, the solution 
processor 40 processes data and optionally can provide the data to the model processor 38. The 

25 model processor 38 can use the data from the solution processor 40 to adjust the one or more 
optical models provided by the model processor 39 to the solution processor 40. 

Also, in alternate embodiments, as indicated by dashed line 37 in FIG. 1, data can be 
shared between the model processor 39 and the data collection data processor 24. The one or 
more optical models provided by the model processor 39 can be used to adjust the electrical 

30 signals provided by data processor 24. This alternate embodiment is described more fully in 
conjunction with FIGS. 5 and 5A. 

It should be appreciated that the same hardware may be used to provide both the data 
collection system 12 and the surface capture system 28. For example, data collection light 
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sources 14 and surface light sources 30 may be provided from the same light hardware which is 
controlled differently so that the light hardware provides light having different patterns (or 
characteristics) depending upon the intended purpose (or function) of the light hardware at a 
particular point in time. For example, if the light hardware were functioning as a data collection 
5 system light source 14, the light hardware would provide light having a first pattern or 

characteristic appropriate for the data collection system function. However, if the same light 
hardware were functioning as a surface capture system light source 30, the light hardware would 
provide light having a second pattern or characteristic appropriate for the surface capture system 
function. The light hardware may be provided, for example, as a programmable light source. 

10 In one particular embodiment, the one or more light sources 14, when used in the data 

collection system 12, generate light in the near infrared (NIR) region and the one or more light 
sensors 20 are adapted to receive the light accordingly. However, in other embodiments, the one 
or more light sources 14 and the one or more light sensors 20 are adapted to transmit and receive 
light having wavelengths above or below the wavelength of NIR light, including visible light and 

15 including infrared light. In one particular embodiment, light sources 14 generate NIR light 
having a wavelength in the range of about 0.630 to 0.950 microns. In another particular 
embodiment, light sources 14 generate IR light having a wavelength in the range of about 0.950 
to 2.0 microns. In another particular embodiment, light sources 14 generate visible light having 
a wavelength in the range of about 0.450 to 0.630 microns. The light provided by the one or 

20 more light sources 14 can be at the same wavelength or a different wavelength than the light 
emitted by fluorochromes described above. One of ordinary skill in the art will appreciate how 
to select a particular wavelength of light (or range of wavelengths) for a particular application. 

In one particular embodiment, the one or more surface light sources 30 generate light in 
the near infrared (NIR), and the one or more surface light sensors 34 are adapted to receive the 

25 NIR light accordingly. However, in other embodiments, the surface light sources 30 and the 
surface light sensors 34 are adapted to transmit and receive light having wavelengths above or 
below the wavelength of NIR light, including visible light and including infrared light. In one 
particular embodiment, light sources 30 generate NIR light having a wavelength in the range of 
about 0.630 to 0.950 microns. In another particular embodiment, light sources 30 generate IR 

30 light having a wavelength in the range of about 0.950 to 2.0 microns. In another particular 
embodiment, light sources 30 generate visible light having a wavelength in the range of about 
0.450 to 0.630 microns. One of ordinary skill in the art will appreciate how to select a particular 
wavelength of light (or range of wavelengths) for a particular application. 
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In another embodiment, the one or more light sources 14, when used in the data 
collection system 12, provide continuous wave (CW) light. However, in other embodiments, the 
one or more light sources 14 are modulated (e.g., in the range of Hz to kHz) or are pulsed (e.g., 
having pulse widths in the range of microseconds to seconds) to enable source multiplexing 
5 and/or background and ambient light separation. Corresponding light detection schemes can also 
be used. 

In other embodiments, frequency domain sources, i.e., intensity modulated light in one or 
multiple frequencies (e.g., in the range of MHz to GHz) or time-domain sources, for example 
pulses of light having different pulse durations (for example, having pulses in the range of 
10 femtoseconds to nanoseconds) and corresponding detection systems can be used. 

In still another embodiment, the one or more light sources 14 provide a planar light 
source, which illuminates at least a portion of the surface of the object 18. In another 
embodiment, the one or more light sources 14 illuminate simultaneously a plurality of spots on a 
surface of the object. One of ordinary skill in the art will recognize that other light patterns can 
15 be also be used. 

It is understood that different projection patterns, such as those described above, may also 
include appropriate masks or spatial attenuation patterns (not shown) to interface to the light 
sensors 16, the light sensors 16 having a dynamic range pertinent to the detection system. With a 
mask, for example, a stray beam of light cannot directly hit and damage or saturate the light 
20 sensors 16. 

Referring now to FIG. 2, a portion 43 of an exemplary tomography system includes one 
or more light sources 46 to provide light 45, which is directed at an object 44. In this particular 
embodiment, the one or more light sources 46 provide a plurality of light beams 45. Also in this 
particular embodiment, the light 45 is directed from the sources 46 toward the object from one 

25 direction. In the exemplary embodiment of FIG. 2 the light sources 46 are shown below the 
object. One or more light sensors 47 are disposed above the object to receive the transmitted 
light 45 having passed through the object 44 and from the object 44 through free space to the one 
or more light sensors 47. It should be appreciated that in some embodiments, the one or more 
light sensors 47, here shown as one light sensor 47, may be moved to a plurality of different 

30 locations as indicated by sensors 47 shown in phantom to receive the light 45. For example, the 
light sources 46 may direct light toward the object from above the object or beside the object in 
which cases the sensors 47 may be stationary or may be moved, scanning the diffuse pattern of 
light propagating inside the object 44 and exiting from the surface of the object 44 at different 

10 
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angles. In other embodiments, a plurality of different light sensors 47 may be disposed in 
particular locations above the object to receive the light 45. It should be appreciated that the one 
or more light sensors 47 can be spaced apart from a surface of the object 44. 

Referring now to FIG. 3, a portion 48 of another exemplary optical tomography system 
5 includes one or more light sensors 47* and one or more light sources 46'. It should be 

appreciated that, in one particular embodiment, a single light source 46* can be moved to project 
light about all surfaces of an object 44', and that the sensors 47' are appropriately positioned to 
receive the light 45 9 having passed through the object 44' and from the object 44' through free 
space to the one or more light sensors 47*. In other embodiments light 45 may be provided from 
1 0 an array of light sources 46* and the light sensors may be provided from an array of light sensors 
47\ 

In still other embodiments, the light sources 46 1 , 46 (FIG. 2) and the light sensors 47', 47 
(FIG 2) are disposed on substantially the same side of the object 44', 44 (FIG. 2). It should be 
apparent that the light sources 46', 46 and the light sensors can be moved or sequenced in tandem 

15 or separately. Also, in yet another embodiment, only the light sensors 47 1 , 47 are moved or 
sequenced while the light sources 46', 46 remain stationary. 

FIGS. 4-5 A are a series of flow diagrams which describe processing performed by a 
system, which may be similar, for example, to the system described above in conjunction with 
FIG. 1, having portions such as those described in conjunction with FIGS. 2 and 3. In FIGS. 4- 

20 5A, rectangular elements are herein denoted "processing blocks" and represent processor 
instructions or groups of instructions (e.g., computer programming code), which may be 
executed by a processing device (e.g., a personal computer, a general purpose computer or any 
other type of suitable processor). Diamond shaped elements, are herein denoted "decision 
blocks," and represent processor instructions or groups of instructions (e.g., computer 

25 programming code) which affect the execution of the instructions represented by the processing 
blocks. 

Alternatively, the processing and decision blocks represent steps performed by 
functionally equivalent circuits such as a digital signal processor circuit or an application specific 
integrated circuit (ASIC). It should be appreciated that the flow diagrams do not depict the 
30 syntax of any particular programming language. Rather, the flow diagrams illustrate the 

functional information one of ordinary skill in the art requires to fabricate circuits or to generate 
computer software or other instruction sets needed to perform the processing required as 
described hereinbelow. It should be noted that many routine program elements, such as 
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initialization of loops and variables and the use of temporary variables are not shown. It will be 
appreciated by those of ordinary skill in the art that unless otherwise indicated herein, the 
particular sequence of steps described is illustrative only and can be varied without departing 
from the spirit of the invention. Thus, unless otherwise stated, the steps described below are 
5 unordered meaning that, when possible, the steps can be performed in any convenient or 
desirable order. 

Turning now to FIG. 4, a flow diagram 49, which illustrates an exemplary process to 
generate a tomographic image with a composite forward model (CFP) method begins as shown 
in processing block 50 in which optical information associated with an object, for example, the 
10 object 18 of FIG. 1, is captured. The optical information can include, for example, information 
provided by the light sensors 20 and the surface light sensors 34 of FIG. 1. 

Processing then proceeds to processing block 52 in which a surface model (e.g., a three- 
dimensional mathematical description) of the surface of the object image is generated. One of 
ordinary skill in the art will understand that there are a variety of ways to optically generate the 
15 three-dimensional mathematical description of the surface. 

First and second optical models are generated as shown in processing block 54. The first 
optical model describes light transmission through the object and the second optical model 
describes light propagation through the surface of the object, in free space about the object and 
characteristics of light sensors used to capture the optical information captured in processing 
20 block 50. 

At processing block 56, a composite forward problem (CFP) is generated as a 
combination of the first and second optical models provided in processing block 54. The CFP 
and optical information associated with the object are used to solve an image problem as shown 
in processing block 58. As described above, the image problem can have the form: 

25 measurements = (theoretical predictions) x (unknown distribution), where the measurements are 
provided from processing block 50 and the theoretical predictions are provided by the first and 
second optical models at processing block 54. Solution of the image problem is further 
described below in conjunction with FIG. 4A. 

Optionally, processing may proceed to processing block 62 in which the first optical 

30 model is adjusted and processing blocks 54-58 are repeated. Once the image problem is solved, 
processing flows to processing block 60 where a tomographic image is generated in displayed. 
Processing then ends. 
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Referring now to FIG. 4A, an exemplary process 63 for generating a tomographic image 
using a composite for a problem (CFP) approach begins by collecting a series of measurements 
as shown in processing blocks 66-74. In particular, calibration measurements, c, are collected in 
processing block 66, background measurements, B, are collected in processing block 68, intrinsic 
5 measurements, I, are collected in processing block 70, fluorescent measurements, F, are collected 
in processing block 72, and surface measurements, R, of the object are collected in processing 
block 74. While the data collected at processing blocks 66, 70, 72, and 74 is collected in the 
presence of light generated at processing block 64, which is directed toward the object, the 
background measurements collected at processing block 68 can be collected in the absence of the 

10 light generated at the processing block 64. 

Calibration measurements can be collected at block 66, where the object, for example the 
object 18 of FIG. 1, can be replaced with an object having known characteristics. For example, 
the object can be replaced with an object having a homogeneous internal structure without 
internal structures. For another example, the object can be replaced with an object having known 

15 internal structures. The calibration measurements thus obtained can be compared later in the 
process 63 with tomographic images generated for the calibration object and the process 63 can 
be adjusted accordingly. The calibration measurements can be collected, for example, with the 
light sensors 20 of FIG. 1. 

Background measurements can be collected in processing block 68 in order to record 

20 ambient light in the absence of the light transmitted at step 64. Background light signals, 

corresponding to the background light can be subtracted or otherwise cancelled from subsequent 
measurements in which with light is transmitted at processing block 64, including from the 
calibration measurements obtained at processing block 66. The background measurements can 
be collected, for example, with the light sensors 20 of FIG. 1. 

25 Intrinsic measurements collected at processing block 70 can include measurements of the 

light that is generated at processing block 64, the light having passed through the object and 
having propagated in free space adjacent to the object 18. The intrinsic measurements can be 
collected, for example, with the light sensors 20 of FIG. 1. 

Fluorescent measurements collected at processing block 72 can include measurements of 

30 fluorescent light generated by fluorochromes from within the object, such as that described in 
conjunction with FIG. 1. The fluorescent measurements can be collected, for example, with the 
light sensors 20 of FIG. 1, with or without the presence of the light transmitted at the processing 
block 64. 
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Surface measurements collected at processing block 74 can include measurements of 
light patterns generated at processing block 64. The surface measurements can be collected, for 
example, with the surface light sensors 34 of FIG. 1. However, as described above in 
conjunction with FIG. 1, in other embodiments, the surface measurements can be collected with 

5 the light sensors 20 of FIG. 1. 

Each of the measurements collected at processing blocks 66-74 are then appropriately 
digitized and filtered as shown in processing blocks 76-84. The data or the information 
measured in processing blocks 66-72 can be used to generate composite measurement vectors, 
M, as shown in processing block 86, which can have the form M=f(m|F, mj, n^B, c), where mi, 

10 m 2 , and m3 are coefficient, any of which can be zero. 

Similarly, the surface measurements of the object having been appropriately digitized and 
filtered at processing block 84 are used to generate a surface model (e.g., a three-dimensional 
(3D) mathematical description) of a surface the object as shown in processing block 88. The 
three-dimensional surface description of the object is then used to provide a first optical model at 

15 processing block 90 as a model of a light field in an object having a shape as described by the 
three-dimensional surface description. The first optical model can assume an object having the 
shape of the actual object to be imaged, the shape provided by the three-dimensional surface 
description, and it can assume an object diffuse to the propagation of light like the actual object. 
However, as described above, the first optical model can assume that the object is homogeneous, 

20 having no internal structures. 

As described above, even for a homogeneous object, light does not propagate in straight 
lines when passing through an object, which is diffuse to the propagation of the light. A variety 
of techniques are known which can provide, as the first optical model, a model of the light field 
in a diffuse object For example, techniques described in PCT Application No. PCT- 

25 US03/17558 filed June 4, 2003, entitled "Imaging Volumes with Arbitrary Geometries in 

Contact and Non-contact Tomography 5 ' can be used. Alternatively, analytical solutions of the 
diffusion equations based on Green's function solutions of homogeneous media, combined with 
first or high-order reflections from the surface elements as when using the Kirchoff 
approximation or the boundary element method can be used. It should be understood that while 

30 the methods described herein represent a first order approximation to the heterogeneous problem, 
the Green f s function solutions can be updated iteratively or, can use a-priori information to 
represent solutions that model heterogeneous media as well. 
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In one embodiment, the first optical model can be calculated based on analytical 
methods, for example the Born or Rytov approximation. However different light field models of 
light in the diffuse medium based on analytical or numerical methods can also be used. 

The three-dimensional surface description of the object provided at the processing block 
5 88 is also used at processing block 92 to provide a second optical model, including, but not 
limited to, a model of light propagation through the boundary of the surface of the object, a 
model of light propagation between the object and light receiver, and a model of characteristics 
of the light receiver. The second optical model is further described in conjunction with FIGS. 
4Aand 6-7E. 

10 The first and second optical models are then combined to generate a composite forward 

problem (CFP) as shown in processing block 94. It should be understood that the CFP can be 
used to predict the output of the light sensors, for example the light sensors 20 of FIG. 1, when 
the object being scanned is homogeneous and has no internal structures. However, when the 
tomographic imaging system 63 is used to scan an object that is not homogeneous or which has 

1 5 internal structures, the measurements collected at the processing block 70 will not agree with the 
prediction made by the CFP. 

The composite forward problem and the composite measurement vectors are then used to 
solve the image problem as shown in block 96. As described above, the image problem can have 
the form: measurements = (theoretical predictions) x (unknown distribution), where the 

20 measurements are provided at processing block 86 and the theoretical predictions are provided as 
the CFP at processing block 94, and where the unknown distribution correspond to internal 
structures in the object. As described above, measurements can also be written as composite 
measurement vectors M=f(miF, m 2 I, m 3 B, c), where mi t m2, m 3 are coefficients, which may take 
zero value. 

25 The first optical model generated at the processing block 90 corresponds to a light field 

description, W. The second optical model generated at the processing block 92 corresponds to an 
operator, T, associated with models of the propagation of light passing through a surface of the 
object, propagation of the light in free space between the object and the light sensors, and 
characteristics of the light sensors. The three-dimensional surface description of the object 

30 generated at the processing block 88 corresponds to a description, S, of the object surface having 
a data set, R. The unknown distribution of optical properties described above, i.e., structures in 
the diffuse object, can be written as a factor, X, corresponding to the unknown distribution 
above. The image problem, therefore, can be written to relate M to X as M = f(miF, m 2 I, m 3 B, c) 
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= q(W(S), T(S), X(S)), or in the more general form as M=g(W'(S), X, S) where q, g are 
appropriate functions that relate the measurements to theoretical predictions and W'(S) is a 
forward descriptor that incorporates light propagation in diffuse and non-diffuse media. The 
functions q or g may be analytical or discrete. In one embodiment described below, M can be 
5 also written as M=T(S)*W(S)*X(S). 

In order to solve the image problem, and in particular to solve for the unknown 
distribution, X, the unknown distribution, X, can be minimized, for example, by an iterative 
solution to the image problem. By another solution method, the function, g, can be inverted to 
extract the unknown distribution, X, from the set of measurements, M. Other methods can also 
10 be used. 

In the CFP method 63, both T and W or W'(s) can be determined theoretically as first and 
second optical models at processing blocks 90 and 92, whereas measurement vectors, M, can 
include measured values provided at the processing clock 86. However, M, in another method, 
referred to as a diffuse surface measurement (DSM) method described in conjunction with FIGS. 
15 5 and 5 A, provides that measurement vectors, M, can be adjusted with theoretical predictions as 
well (i.e., the measurement vectors, M, are scaled by a theoretical calculated function or constant 
associated with an optical model). Similarly, in other embodiments, the operators T , W or 
W'(s), associated with the first and/or the second optical models, can be experimentally 
determined. 

20 Optionally, the model of the light field may be adjusted at processing block 98 and 

processing blocks 90, 92, and 94 may be repeated. Once the image problem is solved at the 
processing block 96, processing flows to processing block 100 where the image is processed 
tomographically, and to processing block 102 where the tomographic image is displayed. 
It should be understood that the composite measurement vectors, M, provided at 

25 processing block 86 are associated with measurements taken with a variety of relative angles 
between light sources and light sensors, in order to provide tomographic image processing at 
processing block 100. 

As described above, in one embodiment, the first optical model generated at block 90 and 
used to calculate the light field, W, in the object can be calculated based on analytical methods, 
30 for example the Born or Rytov approximation. However, in other embodiments, the first optical 
model of the light field in the object can be based upon other analytical or numerical methods. 

To provide the second optical model at block 92, a light source can be considered, at 
position, r, emitting light of wavelength, %u creating an average intensity, U 0 , within an 
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arbitrarily shaped diffuse volume, V, of average absorption coefficient, /* a > and reduced scattering 
coefficient, ft s \ The normalized Born intensity U, (in either absorption/scattering or fluorescence 
mode) is measured by a light sensor at position r d e V , also within the volume. It is given by: 



Uj, U2 are two diffuse photon fields measured experimentally. For example, in 
fluorescence tomography, Ui is the fluorescence field and U2 is the field measured at the 
excitation wavelength. In either absorption tomography or fluorescence tomography, Ui and U2 

10 can be a photon field measured before and after the administration of a contrast agent 

respectively. Or, U2 can be a calibration field measured in a calibration phantom. For different 
cases, the theoretically calculated fields, U3 and U4, can be correspondingly calculated according 
to an assumption that Ui, U2 are the Green's function solutions from the source to the detector, 
or to a voxel, r, respectively. Uo is a background field, which can be turned off or subtracted 

1 5 proportionally to the contact q. This field represents the bleed through signal in fluorescence 
tomography (i.e., U2=Uo and q is the bleed through contact of the filters used). Finally p is an 
"offset" field typically representing some offset in the data due to detection associated offsets or 

CD levels. The photon propagation wave numbers A, =(-v// fll /JD l ) l/2 , ^ =(-v// fl2 /Z> 2 ) 1/2 are 

identical in absorption/scattering mode (kl = k2) and reflect the wave numbers at excitation and 
20 emission wavelengths X\ 9 fa respectively if fluorescence measurements are considered. The 

factors fi aX , fi al , and D x = (3/4 )~\ D 2 = (3/4 y 1 are absorption coefficients and diffusion 

coefficients, respectively. The factor, G t is the system's Green function, which describes light 
propagation from a unit point source to a light sensor and can be written for homogeneous or 
heterogeneous media. The function 0(r) is the unknown quantity reconstructed (this being the 
25 unknown scattering, absorption of fluorochrome mass) and Ao is a multiplicative factor 

associated with system gain factors and light propagation constants such as the speed of light. 
The equation is actually written assuming constant or invariable (known) scattering coefficient 
throughout the medium but can be expanded straightforwardly to include scattering 
inhomogeneities as well. For simplicity, similar light propagation characteristics, e.g., A, « 




(0.1) 



17 



WO 2004/072906 



PCTAJS2004/003229 



and Z>, « D 2 can be assumed, which is a valid assumption when considering the absorption 
scattering problem and is also a valid assumption for fluorescence measurements when there is 
no significant change in tissue optical properties between X\ and Xz 9 which can be true especially 
for dyes in the MR due to the relatively flat absorption spectrum of tissue in this spectral region. 

5 This approximation is done herein for simplification and the methodology works identically if 
different tissue optical properties are considered for X\ and X%. Furthermore, the diffusion 
coefficient, D, is considered to be independent of the absorption coefficient, an approximation 
generally valid in the NIR region. 

While it is computationally very costly to calculate an accurate Green's function for 

1 0 arbitrary boundaries with numerical methods, the Kirchhoff approximation (KA) has been shown 
to be a time-efficient approach. 

In a first assumption, we discretely partition a surface into N plane facets, each one 
having an area, AS b , and a surface normal, n b . We further assume light sensors are at positions 
located outside of the object and in a non-diffusive medium (typically air) where light 

15 propagation can be modeled as straight rays instead of diffuse waves. Therefore, a 

transformation can be determined that describes the contribution of the surface points onto a 
certain light sensor, where light sensors are described as triples, e.g., d = (r d ,n d9 A d ) having 
position, r d , detector normal vector, n d , and aperture, A d . 

In order to generate the above transformation, a variety of optical considerations are 

20 presented in FIGS. 6-7E below. However, first, the diffuse surface measurement (DSM) method 
described above, which also generates the first and second optical models, is further described 
below in conjunction with FIGS. 5 and 5A. 

Referring now to FIG. 5, a flow diagram 120, which illustrates an exemplary process for 
generating a tomographic image with a diffuse surface measurement (DSM) technique begins as 

25 shown at processing block 122 in which optical information associated with an object, for 

example, the object 18 of FIG. 1, is captured. The optical information can include, for example, 
information provided by the light sensors 20 and the surface light sensors 34 of FIG. 1. The 
processing performed at block 122 can be the same as or similar to the processing performed at 
block50 of FIG. 4. 

30 Processing then proceeds to processing block 124, in which a portion of the optical 

information collected at block 122 is adjusted to provide altered measured data referred to herein 
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as diffuse surface measurements (DSMs). The measured data can be altered in conjunction with 
the first optical model as described below. 

Processing then proceeds to processing block 126, in which a surface model (e.g., a three- 
dimensional mathematical description) of the surface of the object image is generated. One of 
5 ordinary skill in the art will understand that there are a variety of ways to optically generate the 
three-dimensional mathematical description of the surface. The processing performed at block 
126 can be the same as or similar to the processing performed at block 52 of FIG. 4. 

The first and second optical models are generated as shown in processing block 128, the 
first optical model to describe light transmission through the object, and the second optical 

10 model to describe light propagation through the surface of the object, in free space about the 
object, and to describe characteristics of light sensors used to capture the optical information at 
block 122. The first optical model can be used to adjust the measured data at block 124 to 
provide the DSM. The first and second optical models can be the same as or similar to the 
optical models generated at block 54 of FIG. 4. 

15 At processing block 130, an image problem is solved. The image problem used in the 

DSM method 120 can have the form: diffuse surface measurements = (theoretical predictions) x 
(unknown distribution), where the altered measurements are provided at processing block 124 
and the theoretical predictions are provided by the first and second optical models at block 28. It 
should be understood that this image problem is very similar to the image problem described in 

20 conjunction with processing block 58 of FIG. 4, and can be solved using the same numerical 

methods. Solution of the image problem is further described above in conjunction with FIG. 4A 
and is not further described here. 

Optionally, processing may proceed to processing block 134 in which the first optical 
model is adjusted in accordance with the image problem and processing blocks 128 and 130 are 

25 repeated. Once the image problem is solved, processing flows to processing block 132 where a 
tomographic image is generated in displayed. Processing then ends. 

Referring now to FIG. 5 A, in which like elements of FIG. 4A are shown having like 
reference designations, a DSM method is shown to have many processing blocks which are the 
same as or similar to processing blocks of the CFP method of FIG. 4A. Here, however, the 

30 processing blocks 152-156 can be different. 

At processing block 152, diffuse surface measurements (DSM) vectors are generated, 
which can have the form M'=f (miF, m2l, n^B, c), where mi, m 2 , and m3 are coefficients. The 
DSM vectors are adjusted versions of the composite measurement vectors, which have the form 
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M=f(miF, m2l, m3B, c) generated at block 86 of FIG. 4A. From analysis above presented in 
conjunction with FIG. 4A, it should be understood that the composite measurement vectors 
represent actual measurements of light collected at blocks 66-72. It should be further understood 
that the composite measurement vectors are used in an image problem at block 96 of FIG. 4A in 
5 order to solve for the unknown distribution, which corresponds to internal structures in the 
object. In the DSM method 150 however, the image problem solved at processing block 156 
receives the diffuse measurement vectors instead of the composite measurement vectors. 

To provide the diffuse surface measurements at processing block 152, the second optical 
model 154 or portions of the second model 154 can provide model information to the processing 

10 performed at processing block 152. The model information provided for this purpose essentially 
moves the measurements collected at one or more of the processing blocks 66-72 by the light 
sensors at a reference position apart from the object to a new reference position on the object, as 
if the measurements had been collected by light sensors on the surface of the object. Therefore, 
the DSM vectors correspond to calculated measurements of the light as if the measurements had 

1 5 been collected by light sensors on the surface of the object. 

At processing block 156 the image problem is solved, having the form: diffuse surface 
measurements = (theoretical predictions) x (unknown distribution). This form has the same form 
as that described above for the CFP method of FIGS. 4 and 4A, and can be solved using the same 
methods. 

20 Referring now to FIG. 6, an imaging system having an aperture, A, corresponding to a 

solid angle, Q 9 generates an image point, R', of a point R lying on a surface, S. Here, the 
surface, S, is a planar surface lying coincident with a focal plane, Zf, associated with the imaging 
system. The image point, R\ lies on an image plane, Z,, and the point, R, lies on the focal plane, 
Zf. Assuming an ideal imaging system, which has a one-to-one correspondence between the 

25 points, R and R\ for those angles that fall within an aperture, A, of the imaging system, 

FCROL =G(R)L, 



where it is assumed thatR f |^ <r> MR\^ , M being a magnification factor. 

The relation between the image at the point, R, that enters the aperture, A, and the image 
30 at the image point, R', is: F (R*) = yF (MR) where y is the overall gain factor (<1). Hereafter, 
the terms R and R' are sometimes used indistinctly, bearing in mind that their relationship is 
R'=MR. Power that goes through the imaging system is represented by: 
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P(*%^=rP{MR)\ m/ (0.2) 

Therefore, we will concentrate on the measurements at the point, R. The "stop," (for 
example, f-stop in conventional photography) of the imaging system is defined by the smallest 
aperture at an entrance to the imaging system. The aperture, A, having an area also denoted as 
5 A, delimits the angle of acceptance of the imaging system . 

Referring now to FIG. 6A, in which like elements of FIG. 6 are shown having like 
reference designations, the plane, 5, having the point R, again is coincident with the focal plane, 
Z& of the imaging system. Therefore, the point, R is in focus in the imaging system (FIG. 6). 
Here, however, the image plane, Zj, of FIG. 6 is not explicitly shown. Instead, the aperture, A, 
10 corresponds to an aperture provided within the imaging system of FIG. 6, and in front of (to the 
left of) the image plane of FIG. 6. 

As described above, the planar surface, 5, is located exactly at the focal plane, Zf. The 
planar surface, S, includes a differential surface area, dB, about the point, R. The differential 
surface area, dB, is shown having a differential surface normal, n. The differential surface area, 
1 5 dB, about the point, R, is associated with a differential surface area, dA 9 at the aperture, A. The 
differential surface area, dA, is shown having a differential surface normal, m. A distance, Za, 
corresponds to a distance from the focal plane, Zf, to the aperture, A, The differential surface 

area, dB, is also shown having a vector, s, in a direction of the differential surface area, dA, and 

* 

the differential surface area, dA, is shown having a vector, u r ».R, in a parallel direction, each 
20 corresponding to an intensity of light passing through the respective differential surface area. 

The aperture, A, can be off-axis with respect to the differential surface normal, n, or it can be on- 
axis. A vector, r\ represents a point corresponding to the point, R, but at the differential surface 
area dA. A vector, O, represents the center of the aperture having area, A. 

A total power measured at the image point, R' (FIG. 6), due to the differential surface 
25 area, dB, is equivalent to the total power, P, at the point, R, (for example, in Watts) irradiated 
into the solid angle, Q (FIG. 6): 

P(R)= J/(r,s)n-sdB^a VseQ (0.3) 
n 

where the vector, n, is the differential surface normal to the differential surface area dB, and 
J(r,s) represents the amount of power that at point r flows within a certain solid angle defined by 
30 unit vector s for each existing wavelength. 7(r,s) is commonly referred to as the specific 
intensity. I(r,s) is shown in FIG. 6A. 

Due to the relationship between r' and the point, R, the differential surface area, dB, may 
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be seen as a "detector area" at the image plane (FIG. 6). In other words, the differential surface 
area, dB, represents an image of the detector area at the focal plane, Zf. The solid angle, ft, (FIG. 
6) corresponds to the aperture, A, (i.e., the entrance pupil, lens, etc.) and to the point, R. In order 
to solve Equation (0.3) it is most convenient to write this solid angle, Q , in terms of the 
5 differential surface area, dA, at the aperture, A, as: 

c/fl= m ' UR -: , flL4 (0.4) 

where m is the differential surface normal at the aperture, A, and where 

u r ,_ R = (r'-R)/|r'-R|, 

10 

where r' is the vector that defines a point in the aperture, A. By using equation (0.4), the total 
power radiated by the differential surface area, dB, is: 

P(R) = f/(R,u r , R )(n • ii^) ( "' B 'g ,) rfftM, (0.5) 
a |r'-R| 

1 5 where the integration is performed over the area of the aperture, A. 
If the aperture, A, has a radius, Ra, then: 

dA = 2nR A dR A and 

20 where Za is a distance from the focal plane (Zf) to the aperture, A. 

Equation (0.5) represents the exact solution within the Radiative Transfer formulation for 
the power, P(fC), (by using Equation (0.2)) measured by a light sensor corresponding to the 
differential surface area, dB\(dB x <r> M 2 dB), where the prime symbol indicates that the 
25 differential surface area dB is the detector area measures at the imaging plane. 

Using the above equations, in conjunction with FIG. 6B below, two case are considered: 
a) when the area of the aperture is small (Z A »R A ), and b) when the area of the aperture is very 
large (Z a «Ra). Below, in conjunction with FIG. 6C, an isotropic light source is considered. 
Referring now to FIG. 6B, in which like elements of FIGS. 6 and 6A are shown having 
30 like reference designations, a small aperture, A, is shown, where the vector, O, represents a 
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center of the aperture, A. Again, the point, R, is in focus in the imaging system (FIG. 6). 

For a small aperture, A, where Z a »Ra, the angular dependence in the integral of 
Equation (0.5) may be approximated to be a constant, yielding: 

P(R) = 7(R,u 0 _ R ) oos0j^j AdB (0.6) 

5 

where O is the vector that represents the center of the aperture, A, and where, 
cos0 = nu o _ R , cosfl^m-Uj^, 

where 9 is the angle between the normal, n, of the differential surface area, dB, and the unit 
vector pointing from O towards R, and 9' is the angle between the surface normal, m, of aperture 

1 0 area, A, and the unit vector pointing from R towards O. 

The above approximation is equivalent, for example, to a light fiber having the aperture, 
A, and surface normal, m, located a 4 distance, O-R, from a radiating source. As an example, the 
case can be considered in which the surface normals, m and n, are parallel (as in the case when a 
planar surface is imaged), and where the differential surface area, dB, radiates as a Lambertian 

1 5 source, /(R,s)=io(R). In this case the power measured at the point, R, is: 

P(R) = / 0 (R)cos 2 9^rdB (0.7) 
r 

where r = |r - 0| is the distance between the points defined by r and O. 

20 For a large aperture, where Za«Ra, most or all light flux that leaves the differential 

surface area dB will be detected. That is, the detected power at the point, R' (FIG. 6), is equal to 
the total power emitted from the differential surface area, dB, and: 

P(R)= J 7(R,s)(n -s)dBdQ = r„(R)dB (0.8) 

where J* is the total flux [e.g., in Watts/cm 2 ] that traverses the differential surface area, dB, 

25 from right to left as shown. 

It should be appreciated that Equation (0.8) is general for any angular dependence of light 
radiated from the image plane, i.e., from the surface, S. For a large aperture, no angular 
dependence introduced due to geometrical considerations is expected in the image plane as long 
as the surface normals, m and n, are parallel. Also, for a large aperture, the in focal plane, Z/, 
30 may be seen as a collection of virtual light sensors (e.g., optical fibers) each having a differential 
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surface area, dB, and each in contact with the surface, S. In the case where the surface normals, 
m and n, are not parallel, the solid angle integral of the above equation will not be in the whole 
hemisphere (2^) + and an angular dependence will appear. 

Referring now to FIG. 6C, in which like elements of FIGS. 6-6B are shown having like 
5 reference designations, again, the point, R, is in focus in the imaging system (FIG. 6). 

In order to understand more clearly the expressions derived above in conjunction with 
FIG. 6B obtained for large and small apertures, a special case is shown and considered where an 
arbitrary aperture, A, is on-axis with the differential surface area, dB. In this case, the total 
detected power is equivalent to: 
10 P(r) = 2^J/(r,s)cos^dBd(cos5) (0.9) 

Q 

where the solid angle, Q, (FIG. 6) is represented by d€l = <tysin 6dQ . The limits of integration 
are cosOe {cos 0^,1} , where cos0 A = Z A /<y]Z A +R] . Assuming an isotropic source, 7(r,s)=/(j, 
Equation (0.9) can be solved to give: 
15 P(r) = xI 0 dB(l-cos 2 9 A ) (0.10) 

and therefore, 

When dealing with the Lambertian approximation, it is shown above that the total 
20 outward flux (i.e., from right to left), J* , is related to the specific intensity as I 0 = J* In . 

Therefore, it should be recognized that when Z a «Ra, equation (0.1 1) is equivalent to equation 
(0.8). Also, since an area, A, of the aperture, A, is A = nR 2 A , equation (0.1 1) can be written as 
P(r) = I 0 AdB/(Zl +R]). In the case where Z a »Ra this reduces to P(r) « I 0 AdBIZ 2 A , thus 

recovering Equation (0.7) for the on-axis case. 
25 Referring now to FIGS. 7 and 7A, in which like elements of FIGS 6-6C are shown having 

like reference designations, the aperture, A, is associated with the imagining system of FIG. 6. 
Here, the focal plane does not lie on a surface, S, of an object being imaged. Therefore, the 
surface, S 9 is substantially out of focus. Also, the in-focus point R is not on the surface, 5, to be 
imaged. 

30 From equation (0.3), the total collected power at the point, R\ is equivalent to the total 
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power radiated by the point, R, into the solid angle, ft , defined by point, R, and the area, A, of 
the aperture, A. 

Referring now to FIG. 7B, in which like elements of FIGS. 6-6C and 7 are shown having 
like reference designations, a differential surface area, dS, which is out of focus, contributes light 

5 energy to point, R\ It will be understood that all of the power that leaves the differential surface 
area, dS, and passes through the point, R, within the solid angle, ft , is detected at the point, R\ 
In order to simplify the calculations, the total power measured at the point, R' (FIG. 7), is 
equivalent to the power that traverses the differential surface area, dB s within the solid angle, Q . 
As described above, this virtual surface represents the image of the light sensor at the point, R\ 

10 The contribution of the differential surface area, dS, can be calculated. It will be 

appreciated that the power radiated must equal the power received. That is, the differential 
surface area, dS, radiates a certain power into the differential solid angle, XI* . From this 
differential radiation, only those values that fall within the solid angle, ft , will contribute to the 
total power. Therefore, the differential power measured at the point, R' (FIG. 7), due to the 

15 differential surface area, dS, may be written as: 

dP(R)= $I f (R,s)m-sdBdn\ Vseft (0.12) 
n 1 

where If is the specific intensity at the focal plane, Zf, the surface normal, m (not shown), is the 
surface normal of the differential surface area, dB, and the solid angle, , is defined by the 
differential surface area, dS 9 and the point, R. 
20 Referring now to FIG. 7C, in which like elements of FIGS. 6-6C and 7-7B are shown 

having like reference designations, another property of the specific intensity, invariance, 
provides that Equation (0.12) above can be written in terms of the specific intensity, J(r,s), at the 
differential surface area, dS: 

rfP(R)= J/(r,s)n'SdStfft', Vseft, (0.13) 

25 where the solid angle, <5ft 1 , is defined by the differential surface area, dB, and the point r is the 
surface point at which the differential surface area, dS, is located. It should be recognized that 
the only values of, s, that contribute to the solid angle integral in Equation (0,13) are those that 
fall within the solid angle, ft (FIG. 7B). It should be appreciated that the above equation 
corresponds to the total power received at the point, R' (FIG. 7), due to a differential surface 

30 area, dS. 

Referring now to FIG. 7D, in which like elements of FIGS. 6-6C and 7-7C are shown 
having like reference designations the total power takes into account the complete surface, for 
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those values that fall within the solid angle, Q (FIGS. 7 and 7B) through the differential surface 
area, dB (FIG. 7B). This is schematically represented in Fig. 7D. Hie total power received at 
the point, R 9 (FIG. 7), is therefore: 

P(R) = JdP(R) = \dS J/(r,s)n -sdQ.\ Vs e Q (0.14) 
s s n» 

5 Referring now to FIG. 7E, in which like elements of FIGS. 6-6C and 7-7D are shown 

having like reference designations, for a small aperture, where Za»Ra, and an out of focus 
surface, 5, it is convenient to return to the contribution of a single differential surface area, dS, 
given in Equation (0.13). Since the solid angle in this case is a delta function, we may rewrite 
Equation (0.13) as: 

10 

ctf > (R) = /(r,s)n•s^(s-u r „ 0 )J(s-u 0 ^ R )^Q , , (0.15) 
where the delta functions ensure that the direction of energy that passing through the points, R 
and O, is taken into account. Rewriting the solid angle, dCl\ in terms of the differential surface 
area, dB, we obtain: 

15 dP(R) = /(r,s)n.sJ(s~u r ^ 0 )J(s~u 0 ^J^^^, (0.16) 

which is equivalent to: 

dP(R) = /(r,u r _ R )*(u r _ R -u^Xu^ -n)(u r _ R ^) ^fp (0-17) 

Introducing this expression into Equation (0.14)we obtain the total power received at the 
point, R J , (FIG. 7) is: 

20 P(JR) = /(r RO ,u M )[u RM) n(r R0 )][u 0 _ R -m] ^ (0.18) 

l r RO~ R l 

where tro is the point where the surface, 5, and the line given by RO intersect. That is, from 
Equation (0. 1 8) it is clear that in the case of a small aperture only one element of the surface, S, 
will contribute to the intensity at a certain pixel. What this represents is that for a very small 
aperture, A, all of surface, S, will be in focus. The differential surface area, dB, represents the 
25 image of the light sensor at the focal plane, Z/ 

For a large aperture, where Ra»Za, all surface points on the surface, 5, contribute to the 
power measured at the point, R' (FIG. 7). Rewriting the solid angle in terms of the differential 
surface area, dB. 
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P(R) = J/(r,u r . R )[n u r _ R ][ro j— ^ (°- 19 ) 

That is, for all points out of focus, there is a complete de-blurring. Therefore, very large 
apertures have a very sharp focal plane, Zf, or focus depth. 

It should be recognized that for all of the expressions where I(r,s) has not been defined 
5 explicitly, the above expression includes all high-order angular components of the light at the 
point, r (FIG. 7D). For an arbitrary surface, S, being imaged we may assume the transformation 
of a general point, rt,, (equivalent to point r) within the surface, S, and a light detector at a 
general point, r d (equivalent to point R) within the focal plane Zf as: 

r(r„r,) = £(r 6 ,rj/(n,m)^, (0.20) 

10 

In Equation (0.20) n and m represent the surface normal, n, and the detector normal, m, 
respectively, as before. Function £ is a visibility factor that discards surface points not visible 
from the light sensor, whereas function /includes the angles of acceptance of the aperture, A, 
with regards to each surface component, defined by n and m. Function T(r^r rf ) includes all the 
1 5 high-order angular contributions of light and may therefore be written as 

r(r 6 ,rj = Yja n r in) (r by r d ) , where T in) (r b ,r d ) is the n-order angular contribution, i.e., a function 

n 

of the n-order Legendre Polynomial, and a n is the contribution of this order. 

Using Equation (0.20) we may relate a set of light intensities Ufart) at the diffuse/non- 
diffuse air-tissue interface delineated by the surface, S , to non-contact light intensity 
20 measurements U nc (r 5 , r d ) obtained from a free-space light sensor at r d £ V , i.e., 

^(r.,r < ,) = 2 r (r rf ,r i )C/(r J ,r,)A5 (0.21) 

r h eS 

where T(r b ,r d ) is a matrix. In the case where there is a one to one correspondence (surface in 
focus), T(r 6 , r rf ) will be a diagonal matrix. Equation (0,21) may also be envisaged as a 
convolution of the surface values with the transformation function r(r b ,r d ) . 
25 The transformation of Equation (0.21) makes no assumptions regarding the position, 

orientation, or area of the surface elements or the light detectors but explicitly models these 
parameters through Equation(0.20). 

As we consider a finite set of detectors and surface elements, Equation (0.21) can be 

rewritten as a matrix equation where V m e tf |dct|H , T e R M ® and U e # H describe an 
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imaging system with \det\ detectors, \src\ sources and \S\ surface elements. Each matrix 
element (b,s) in U contains the average intensity on a surface element r b e S created by a 
source r, as described by Equation (OA). Transforming the volume integral into a sum of 
volume elements (voxels) of volume AV r leads to: 



C/ 0 (r„r,)r(r A ,rJ 
where the sum is performed over all the elements within V. 

Equation (0.22) can be considered to be the product of a matrix W e with a vector 
c e Jl M , where the matrix contains all the factors in Equation (022) except for c. Thus, the 
10 complete discrete system can be written U wc = Wc . As W tends to be very large, the inversion 
cannot be performed directly. Instead, iterative techniques, such as the algebraic reconstruction 
technique with randomized projection order (R-ART) can be used, conjugate minimization 
methods, or other methods. 

All references cited herein are hereby incorporated herein by reference in their entirety. 
1 5 Having described preferred embodiments of the invention, it will now become apparent 

to one of ordinary skill in the art that other embodiments incorporating their concepts may be 
used. It is felt therefore that these embodiments should not be limited to disclosed embodiments, 
but rather should be limited only by the spirit and scope of the appended claims. 

20 What is claimed is: 



28 



WO 2004/072906 



PCT/US2004/003229 



CLAIMS 

11. An system for optical tomography of an object, comprising: 

2 a data collection system to collect light having traveled through the object and in free 

3 space about the object and to generate one or more signals in response to the light; 

4 a surface-capture system to receive surface light reflected from a surface of the object and 

5 to generate a surface model of at least a portion of the surface of the object in response to the 

6 surface light; 

7 a model processor coupled to the surface-capture system to provide one or more optical 

8 models associated with the surface model; and 

9 a solution processor coupled to the data collection system and to the model processor to 

10 provide a tomographic image of the object in response to the one or more signals and to the one 

1 1 or more optical models. 

1 2. The system of Claim 1 , wherein the data collection system comprises: 

2 one or more light sources for emitting light directed to propagate through the object, the 

3 object diffuse to propagation of the light through the object; 

4 one or more light sensors disposed to receive the light which has passed through the 

5 object, wherein at least one of the one or more light sources and the one or more light sensors are 

6 disposed apart from the object; and 

7 a data processor to generate the one or more electrical signals in response to the light 

8 received by the one or more light sensors. 

1 3. The system of Claim 2, wherein the one or more light sensors corresponds to one light' 

2 sensor. 

1 4. The system of Claim 2, wherein the one or more light sources correspond to one light 

2 source. 

15. The system of Claim 2, wherein the one or more light sensors are adapted to receive 

2 fluorescent light emitted from fluorochromes within the object. 

1 6. The system of Claim 2, wherein the one or more optical models include a first optical 
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2 model to provide a description of light in the object assuming that the object is homogeneous and 

3 diffuse to light propagation, and a second optical model formed by modeling at least one of: 

4 light propagation associated with the light passing through the surface of the object, propagation 

5 of the light in free space between the object and the one or more light sensors, and the one or 

6 more light sensors. 

17. The system of Claim 6, wherein the model processor comprises a composite-forward- 
2 problem (CFP) processor for combining the first optical model and the second optical model. 

18. The system of Claim 7, wherein the solution processor provides an image solution by 
2 relating the CFP to the one or more electrical signals. 

1 9. The system of Claim 6, wherein the data processor comprises a diffuse surface 

2 measurement (DSM) processor for adjusting the one or more electrical signals in response to the 

3 second optical model output. 



1 1 0. The system of Claim 2, wherein the surface-capture system comprises: 

2 a surface light source for emitting surface light directed toward the surface of the object; 

3 the one or more light sensors disposed to receive the surface light reflected from the 

4 surface of the object to provide one or more image signals; and 

5 a surface-capture processor to generate the surface model in response to the one or more 

6 image signals. 

1 11. The system of Claim 1 , wherein the surface-capture system comprises: 

2 a surface light source for emitting surface light directed toward a surface of the object; 

3 a surface light sensor disposed to receive the surface light reflected from the surface of 

4 the object to provide one or more image signals; and 

5 a surface-capture processor to generate the surface model in response to the one or more 

6 image signals. 



1 12. The system of Claim 1, wherein the light is in near infra-red (NIR) light having a 
1 wavelength between about 0.630 and 0.950 microns. 
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1 13. The system of Claim 1, wherein the light is infra-red (IR) light having a wavelength 

2 between about 0.950 and 2.0 microns. 

1 14. The system of Claim 1, wherein the light is visible light having a wavelength between 

2 about 0.450 and 0.630 microns. 

1 15. The system of Claim 1 , wherein the solution processor is adapted to provide the 

2 tomographic image for the object having an arbitrary shape. 

1 16. The system of Claim 1, further comprising: 

2 an image processor coupled to the solution processor to process the tomographic image; 

3 and 

4 a display coupled to the image processor to display the tomographic image. 

1 17. A method of imaging an object with optical tomography comprising: 

2 capturing light having traveled through the object and in free space about the object; 

3 generating one or more signals in response to the light; 

4 generating a surface model of at least a portion of a surface of the object; 

5 generating one or more optical models associated with the surface model; and 

6 generating a tomographic image of the object in response to the one or more signals and 

7 to the one or more optical models. 

1 18. The method of Claim 17, wherein the capturing light comprises: 

2 emitting light with one or more light sources, the light directed to propagate through the 

3 object, the object diffuse to propagation of the light through the object; 

4 receiving the light with one or more light sensors, the light having passed through the 

5 object, wherein at least one of the one or more light sources and the one or more light sensors are 

6 disposed apart from the object; and 

7 processing outputs of the one or more light sensors to generate one or more electrical 

8 signals in response to the light. 

1 1 9. The method of Claim 1 8, wherein the one or more light sensors correspond to one light 

2 sensor. 
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1 20. The method of Claim 1 8, wherein the one or more light sources correspond to one light 

2 source. 

1 21. The method of Claim 1 8, further comprising receiving, with the one or more light 

2 sensors, fluorescent light emitted from fluorochromes within the object. 

1 22. The method of Claim 1 8, wherein generating the one or more optical models comprises: 

2 modeling, as a first optical model, the light in the object assuming that the object is 

3 homogeneous and diffuse to light propagation; and 

4 modeling, as a second optical model, at least one of: light propagation associated with 

5 the light passing through the surface of the object, propagation of the light in free space between 

6 the object and the one or more light sensors, and the one or more light sensors. 

1 23 . The method of Claim 22 further comprising combining the first optical model and the 

2 second optical model to provide a composite-forward-problem (CFP). 

1 24. The method of Claim 23, wherein generating the tomographic image comprises providing 

2 an image solution by relating the CFP to the one or more electrical signals. 

1 25 . The method of Claim 22, further comprises adjusting the one or more electrical signals in 

2 response to the second optical model output to provide a diffuse surface measurement (DSM). 



1 26. The method of Claim 1 8, wherein generating the surface model comprises: 

2 emitting surface light with a surface light source, the surface light directed toward the 

3 surface of the obj ect; 

4 receiving the surface light reflected from the surface of the object with the one or more 

5 light sensors to generate one or more image signals; 

6 generating the surface model of the at least a portion of the surface of the object in 

7 response to the one or more mage signals. 

1 27. The method of Claim 17, wherein generating the surface model comprises: 

2 emitting surface light with a surface light source, the surface light directed toward a 
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3 surface of the object; 

4 receiving the surface light reflected from the surface of the object with a three- 

5 dimensional light sensor to generate one or more image signals; 

6 generating the surface model of the at least a portion of the surface of the object in 

7 response to the one or more image signals. 

1 28. The method of Claim 1 7, wherein the light is near infra-red (NIR) light having a 

2 wavelength between about 0.630 and 0.950 microns. 

1 29. The method of Claim 17, wherein the light is infra-red (IR) light having a wavelength 

2 between about 0.950 and 2.0 microns. 

1 30. The method of Claim 17, wherein the light is visible light having a wavelength between 

2 about 0.450 and 0.630 microns. 

1 31. The method of Claim 1 7, wherein the generating the tomographic image of the object 

2 provides the tomographic image for the object having an arbitrary shape. 
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